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LATERAL STABILITY CHARACTERISTICS AT LOW LIFT BETWEEN 
MACH NUMBERS OF O .85 AND I.I 5 OF A ROCKET-PROPELLED 
M3DEL OF A SUPERSONIC AIRPLANE CONFIGURATION 
HAVING A TAPERED WING WITH CIRCULAR-ARC 
SECTIONS AND kO° SWEEPBACK 
By Charles T. D'Alutolo and Allen B. Henning 

SUMMARY 


A rocket -propelled model of a supersonic airplane configuration 
eii 5 )loying a 40° swepthack wing having circular -arc sections was flight 
tested in the Mach number range of O .85 and I.I 5 to obtain lateral sta- 
bility data at low lift from the lateral response to an inrpulsive rudder 
deflection and to evaluate the test and analysis techniq.ue. 

Time histories of the Dutch-roll oscillations indicate that the 
model was statically stable throijghout the Mach number range but was 
dynamically unstable over a small range of Mach numbers. A low -lift 
bviffet was experienced below a Mach number of 0.95» 

The time-vector method applied to the recorded Dutch-roll transient 
oscillations provided a usefvil method for the determination of the lateral 
stability derivatives. The results as obtained from the vector analysis 
indicate that the directional stability and effective dihedral increased 
with increasing Mach number, the danping -in-yaw was low and over a small 
region of Mach number was unstable, and the rate of change of rolling- 
moment coefficient with yawing-angular -velocity factor was negative at a 
Mach number of O .89 and increased positively to a large positive value 
at a Mach number of 1.1. The lateral -force derivative increased with 
increasing Mach number until a Mach number of O .98 and then decreased to 
the limit Mach number of the test. 

Conparlsons between the rocket-propelled-model test data, wind-tunnel 
data, and estimates were made in order to evaluate the test technique. 
These comparisons indicated that the time-vector method allows the 
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determination of the static lateral stability derivatives to the same 
order of accuracy as does wind-tunnel techniq_ues . 


INTRODUCTION 


The Pilotless Aircraft Research Division is conducting a flight 
investigation to determine the longitudinal and lateral stability ch£.rac— 
teristics at low lift in the transonic speed range of a supersonic air- 
plane configuration having a tapered wing with circular -arc sections and 
j^o sweepback. The longitudinal stability characteristics of the coiifig- 
uration are presented in reference 1 and the present paper contains ihe 
results from a flight to determine the lateral stability characterisi.ics. 
The Mach number range covered in the present test was from O .85 to I.I 5 

and corresponds to a Reynolds number range of 7*5 X 10^ to 9*8 X 10^ j 
respectively. The model was flown at the Langley Pilotless Aircraft 
Research Station at Wallops Island, Va. 

Stability derivatives were determined by application of the time^- 
vector method (see refs. 2 to 6) and these derivatives are con^jared i-lth 
other tests and with the estimated values of the stability derivatives 
in order to evaluate the test techniq.ue. The physical motions experj.- 
enced by the model may or may not be the same as those experienced bj- 
the airplane since the model and airplane would have different mass £Jid 
inertia characteristics . 


SYMBOLS AND COEFFICIENTS 


Throiaghout this paper the forces and moments acting on the mode], are 
referred to the body system of axes, which are defined as an orthogonal 
system of axes Intersecting at the airplane center of gravity, in which 
the Z-axis is in the plane of symmetry and perpendicular to the X-axis. 

The X-axis is in the plane of symmetry, and the Y-axis is perpendicu].ar 
to the plane of symmetry. A diagram of these axes showing the positive 
direction of forces, moments, and angles is presented in flgiire 1. 55ince 
aerodynamic derivatives are usually available relative to the stabllf.ty 
system of axis, a diagram showing the stability system of axis is inc:luded 
in figure 1 for reference purposes. 

Reference 7 gives expressions that can be losed to transpose the aero- 
dynamic derivatives from one system of axes to the other . Angular rela- 
tionships in flight for the stability and body systeit^ of axes indiciiting 
positive directions of angles is presented in figure 2. 
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The symbols and coefficients are defined as follows; 


A 

HU 

H/s 

^T/g 


a 


b 


c 


'W 


aspect ratio 

acceleration along X reference axis as obtained from accel- 
erometer, positive forward 

acceleration normal to X reference axis as obtained from 
accelerometer, positive \jp 

acceleration along Y reference axis as obtained from accel- 
erometer, positive to the right 

total damp ing factor (logarithmic decrement of Dutch-roll 
oscillation defined as being a positive number for a damped 
oscillation) 

wing span, ft 

wing chord, ft 

mean aerodynamic chord of wing, ft 


c^ mean aerodynamic chord of vertical tail, ft 

D differential operator, d/dt 

2 

I^ moment of Inertia about body X-axis, slug -ft 

I2 moment of Inertia about body Z-axis, slvig-ft^ 

Ij^ product of inertia referred to body axis (positive when the 

positive direction of the X principal axis is inclined 
below reference axis, that is, when e is positive) 

\\^Z - Ix]‘'^a.n2e 

L rolling moment, ft-lb 

M' pitching moment, ft-lb 

N yawing moment, ft-lb 

M Mach number 


m mass of model, slugs 
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period of Dutch-roll oscillations, sec 

rolling angular velocity about X-axis , radians/ sec 

total stagnation pressure, Ib/sq. ft 

dynamic pressvtre, ipV^, Ib/sq ft 

Reynolds number 

yawing angular velocity about Z-axis, radians/sec 
total wing area, sq. ft 

one -half thickness of airfoil at aileron hinge line, as shown 
in figure 3(12) • 

time, sec 

velocity, ft/sec 

model weight, lb 

air density, slugs/cu ft 

relative density factor, -S_ 

pSb 

undamped natural circular freq.uency, (o)2 + 
radians /sec 

frequency of the Dutch-roll oscillation, radians/sec 

angle in plane of symmetry, measiired from projection of rela- 
tive wind to fuselage reference axis, deg or radians (see 
fig. 2) 

angle of sideslip, measured from relative wind to fuselage 
reference axis, deg or radians (see fig. l) 

angle between principal longitudinal axis of inertia and the 
longitudinal body axis, deg, (see fig. 2.) 

flight path angle, deg, (see fig. 2) 

angle between reference axis and horizontal axis, deg, (see 
fig. 2) 
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C 


inclination of principal longitudinal axis of inertia with 
respect to flight path, deg, (see fig. 2) 

angle of roll, radians 

angle of yaw, radians 

rudder deflection, deg 


phase angle, deg 
trim lift coefficient 

lateral -force coefficient. Lateral force/q.S 
yawing -moment coefficient. Yawing moment/qSb 
rolling -mcanent coefficient. Rolling moment/qSb 
directional stability derivative, — =, per radian 

ap 

lateral -force derivative, - — y per radian 


effective dihedral derivative. 


^ — , per radian 

^P 


rate of change of yawing-moment coefficient with rolling- 
angular -velocity factor, — per radian 

2V 


rate of change of lateral -force coefficient with rolling- 


angular -velocity factor. 





per radian 


^7 

dajnping-in-roll derivative, — per radian 

2W 
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rate of change of yawing -moment coefficient with yawing- 
angular -velocity factor, — — , per radian 

sik 

2V 

rate of change of lateral-force coefficient with yawing- 

^Cy 

angular-velocity factor, — -, per radian 

2V 

rate of change of rolling-moment coefficient with yawing- 

bC. 

angiilar -velocity factor, — radian 

S— 

2V 

rate of change of yawing-moment coefficient with rate of 

dc 

change of angle-of -sideslip factor, — per radian 

skk 

2V 




'la 


rate of change of lateral -force coefficient with rate of 

SCy 

change of angle-of -sideslip factor, — —f per radian 

2V 

rate of change of rolling-moment coefficient with rate cf 

SCj 

change of angle -of -side slip factor, ; , per radian 

dE 


The symbol [j| represents the absolute magnitude of j and is 
always taken to he positive. A dot over a variable indicates the first 
derivative of the variable with respect to time. Two dots indicate the 
second derivative. Phase angles are indicated by subscript notation 
as <5 as the phase angle between sideslip and rolling angular ac:cel- 

eration. The second subscript symbol is ixsed as the reference. A posi- 
tive sign associated with the phase angle indicates that the first sub- 
script symbol leads the reference, whereas a negative sign indicates 
that the first subscript symbol lags the reference. 
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MODEL AND TEST DESCRIPTION 
Model Description 


The general arrangement of the model ^ details of wing and tall, and 
details of rudder are shown in figure 5, and the geometric and mass char 
acteristics of the model are given in table I. Photographs of the model 
and the model-booster combination are shown in figure k. 

The model fuselage was a body of revolution of fineness ratio 
consisting of a cylindrical center section, ogival nose and tail sec- 
tions, dorsal fin, canopies, and fuselage skid. Construction of the 
fuselage was principally of aluminum covered with magnesium skin. 

The nose section contained the telemeter; the center section con- 
tained the power section and wing mount; and the tail section contained 
the rudder -pulsing mechanism. 

The wing of the model was made of steel and had 10 -percent circular 
arc airfoil sections perpendicular to the quarter -chord line and incor- 
porated k0° of sweepback at the quarter -chord line with 5° positive 
dihedral. The wing was set at 5° incidence with respect to the fuselage 
reference line and was modified to simulate slab-sided ailerons having a 
trailing-edge thickness of one-half the aileron thickness at the hinge 
line. The ailerons were of 25-percent span and were constructed (or 
set) at 0° deflection. Tip fairings were placed on the wing in order to 
house some of the Instrumentation. 

The horizontal tail was similar to the wing in plan form but had 
MCA 65-008 airfoil sections parallel to the fuselage reference line and 
was constructed of aluminum. An incidence of 2° trailing edge down 
relative to the f\iselage reference line was set in the horizontal tail 
so that the model wo\ald have reasonable trim values. (See ref. l) . 

The vertical tail had MCA 27-010 airfoil sections parallel to 
fuselage reference line at the root and MCA 27-008 airfoil sections 
parallel to fuselage reference line at the tip and was constructed of 
aluminum. An aluminm rudder, the details of which are shown in fig- 
ure 5 (c), was incorporated on the vertical tail and was used to disturb 
the model in yaw. 

The model rudder -pulsing mechanism was designed to deflect the 
rudder impulsively between 0^ and 25^ by means of a hydraulic servo 
system and then allow a dwell time before the rudder was deflected again 
The time required to move the rudder from 0° to the maximimi deflection 
of 25° and then back to 0° was 0.07 second and the dwell time was about 
0.54 second; thus, a con^lete cycle occurred every 0.6I second. 
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The model was boosted to a Mach nvuiiber of l.l 8 by an external 
A£L Deacon rocket motor* Upon burnout of this rocket motor ^ the model 
separated from the booster and coasted through the test speed range 
(data obtained from M = I.I5). The model-booster combination was 
launched from a mobile launcher at an angle of about as shown ki 
figure 4 (c) , 


Instrmentation 

The model contained a standard twelve -channel MCA telemeter. 
Measurements were made of the normal, longitudinal, and transverse accel- 
erations near the center of gravity of the model and the normal and 
transverse accelerations in the nose of the model. Normal accelerabions 
of each wing tip were measured in order to determine the rolling angular 
acceleration of the model. Rolling angular accelerations were also 
measured by an angular -accelerometer -type instriiment. The angle of attack 
q nti angle of sideslip were measured by a vane-type instrument located on 
a sting forward of the nose of the model, whereas total pressure was 
measured by a tube located on a small strut mounted on the underside of 
the model near the nose and rudder deflections were measured by a control- 
position pickup. 

The position of the model in space was determined by use of a modi- 
fied SCR 584 tracking radar set and the velocity of the model was obtained 
by use of the CW Doppler velocimeter radar set. Atmospheric data were 
obtained from a radiosonde released Immediately before the model flight. 
Fixed and tracking motion-picture cameras were used to observe the con- 
dition of the model during most of the flight. 


Preflight Test 

Prior to flight testing, the model was suspended by shock cords and 
shaken by of an electromagnetic shaker. This preflight test was 

performed in anticipation of the model ejcper fencing hlgh-freq.uency oscil- 
lations during the flighty thereby, an explanation of these high— freq^uency 
oscillations was possible. The frequencies recorded in the shake tests 
of the model together with the approximate nodal lines are tabulated on 
the following page. 
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Mode 

Nodal line 

Frequency, 

cps 

Wing first 
bending 

Vane / j 


Horizontal- 
tail first 
bending 

/] 

86 

Wing second 
bending 

n 

Vane / / 

176 

Horizontal- 
tail second 
bending 

/7 

356 

Wing torsion 

/7 

Vane / ^ / 

262 


Flight Test 

The model was boosted to a Mach mimber of I.18 and upon b^I^nout of 
the booster rocket motor the model separated from the booster. Dviring 
the boosted phase of the flight, the rudder -pulsing mechanism was 
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inoperative and was not allowed to operate until the model was completely 
separated from the booster (M = I.I 5 ) . After completely separating from 
the booster, the model was disturbed in yaw by periodic pulses of the 
rudder and time histories of the resulting model motions were obtained 
by means of the NACA telemeter and instrument system. 

The flight conditions of the model are presented in figure 5 where 
the variation of air density, velocity, dynamic pressixre, and relati-ve 
density factor with Mach number are shown. These q.uantities are pre- 
sented so that a possible correlation of the data obtained from this 
test with data obtained from other tests may be made. The range of ihe 
Reynolds niimbers of the present test is shown in figure 6. 


ACCURACY AND CORRECTIONS 
Accuracy 


The estimated probable errors in the basic quantities measured are 
shown in table II. The lateral stability derivatives Cv , C„ , C 7 , 

P 3 P 

Cl , and C„ - C„, are dependent upon some of or all these measured. 

‘r “r 

quantities. The probable error in any of the above derivatives due to 
all the probable errors in table II was determined by the method shovn 
in reference 2 and it is felt that the accuracy of the derivatives 
reported in this paper are of the same order as those reported in ref- 
erence 2. That is, at M = 1.1 and M = 0.9 the accuracies are, respec- 


tively: 5 and 5 percent for Cy 7 6 and 8 percent for Cj^ and Cj 

P P 

13 and 16 percent for C^ , and I5 and 25 percent for Cj^ - Cj^«. 

r ip 


p' 


It is believed that the data presented in this report provide a good 
indication of the variation of the stability derivatives with Mach number, 
flurl the absolute values of these derivatives are at least as accuirate or 
better than indicated above. 


Corrections 

Since it was impossible to mount the accelerometers exactly at ihe 
center of gravity of the model, the accelerometer readings had to be 
corrected so that the accelerations of the center of gravity of the n.odel 
could be determined. These corrections consisted of the linear and 
angular acceleration effects as well as angular velocity effects on ihe 
accelerometer readings due to the model motions. 
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The angles of attack and sideslip as measured in front of the model 
were corrected to the model center of gravity by considering flight path 
curvature effects, see reference 8. 

Frequency -response corrections to all instr\jments except the angular 
accelerometer were not necessary since the model natirral frequency was 
less than 4 percent of the instrument natural frequency. Frequency- 
response corrections, however, were applied to the angular -accelerometer 
readings since the model natural frequency was as high as 20 percent of 
the angular -accelerometer natural frequency. These corrections were 
made by lise of standard frequency -response charts and resulted in cor- 
rections to the phase angle only. The frequency-response corrections to 
the phase angle between the rolling angular acceleration and sideslip 
angle amoimted to 9-6° at M = 0.89 and 16.5° at M = 1.10. 


AMLYSIS 


The method used in this paper to analyze the data and to determine 


and C„ - C. 


the lateral stability derivatives ^n^ “ ^n^ 

based on the concept of rotating vectors. This concept as applied to 
airplane dynamics was first formulated by Mueller (see ref. 3) aJ^d 
extended by others (see refs. 4 to 6) and may be briefly described as 
follows ; 


IS 


A fundamental property of rotating or time vectors is the relation- 
ship between the vector and its derivatives or its integrals at a given 
instant of time. This relationship may be illustrated by considering 
the oscillatory motions resulting from a yaw disturbance which are 
asstimed to be given by 


3 



cos cut 


( 1 ) 


9* - |p[ Vp) 


Differentiation of equation (5) gives 


^ = Ja.^ + oi^ 3o®"^^ cos^cot + + 90° + tan"^ ^ (4) 
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whereas the second differentiation of equation (3) yields 

<f) = (a^ + cos^ojt + + 180° + 2 tan (5) 

where o) is the frequency at which these vectors are rotating and Pq 
is the initial value of the sideslip angle. 

Differentiation of equations (l) and (2) yields similar expressions. 

Coii5)arisons of equations (3) and ( 4 ) show that the an5>littide of the 
first derivative of a time vector is equal to the amplitude of the time 
vector at the same instant of time multiplied by the mdamped natural 
circular frequency of the vector^ and that the phase of the first deriv- 
ative leads the time vector by 90° + tan"^ ^ where tan"^ ~ is referred 

to as the damping angle. 

Coitparison of eqijations (3) and (5) shows that the amplitude of the 
second derivative of a time vector is equal to the amplitude of the time 
vector multiplied by the square of the mdamped natural circular fre- 
quency of the vector, and that the phase of the second derivative leads 
the time vector by 180° plus twice the damping angle. 

In a like manner, the amplitude of the integral of a time vector 
may be shown to be equal to the time vector divided by the \jndamped 
natural circular frequency, whereas the phase of the integral of a time 
vector lags the time vector 90° pl\is the damping angle. 

Another important property of the concept of time vectors is the 
requirement that the vector polygon representing any degree of freedom 
of a system must close; thereby the determination of only two unknown 
quantities in the degree of freedom is allowed. This property will be 
explained in detail imder the discussion of the solution of the lateral 
equations of motion. 

The complete lateral equations of motion about fixed body axis (see 
fig. 1) are: 

Lateral force: 

^(Dp + Di|a - oD^) = Cy^P + ^ ^ ^ + 


— (rf cos 9 + \|r sin 0) 
qS 


( 6 ) 
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Rolling moment: 

d2o( - D2,|r 

qSb qSb ^ 


= C 


1 3 + 

ip 


-L + 

ip 2V ^ 


4- Dik + C 
ir 2V ^ 


ip h 


(T) 


Yawing moment: 

^ ^ S I? * ""r I? '’♦ * S If 


( 8 ) 


Assumptions were made concerning some of the terms in these equa- 
tions in order to facilitate analysis of the data. These assumptions 
are: In the lateral -force equation all the aerodynamic terms are com- 

bined into one term referred to as Cy or the total lateral force. 

This assumption is valid since the total lateral force was measured by 
a transverse accelerometer and includes the contributions of rolling 
angular velocity, yawing angular velocity, and sideslip angle. It was 


further assimaed that C 


Y was eqixal to ' 


The validity of this 


assumption is shown from the time history of the model motions where Cy 

is within the limits of accuracy of the test l 80 ° out of phase with the 
sideslip angle p. It was also assumed that the gravity terms in the 
lateral -force eqmtion may be neglected. This assumption is valid for 
rocket-propelled models since the models have low wing loadings and are 
flown through rather dense air at high speeds so that the values of the 
gravity terms are very small. In the rolling -moment eqmtion it was 
assmed that C^^ = 0 , whereas in the yawing -moment equation it was 

assimied that \|r = -P, so that C^j^ and Cj^, may be combined into one 

term referred to as the damp Ing-in -yaw derivative. 

By considering the above assimptions , the lateral equations of 
motion sire now written in the form used to analyze the data by the time- 
vector method. 


Lateral force: 


mv /lpl l\jf| 1^1 IpI 

asl|3l |3l ■ |3|r 131 IPI " 


( 9 ) 


Rolling moment: 


% \'i\ ^xz iVl 

qSb |p| qSb |Pl 


-C7 i-JiL = o (10) 

h |3| h 2V |P| V 2V |P| 


CONFIDENTIAL 


Ik 


CONFIDENTIAL 


NACA EM L55A51 


Yawing moment; 

^ ^ M . M . i M -fc_ - C„.\ ii|- = 0 (11) 

<lSb Ipl qSb |3( ""p IpI ^ 2V |3| py 2V IpI 

Notice that each of the equations was divided by p which is considered 
as the unit vector in the time-vector plots and the various amplitude 
ratios were determined from the data. 

From preflight measurements, the mass, geometric, and inertia char- 
acteristics of the model were determined; from the radar, radiosonde , and 
telemeter -instrument system, the velocity, dynamic pressure, Mach number, 
lateral force, angle of attack, angle of sideslip, rolling angular accel- 
eration, and yawing angular acceleration were obtained. Thus, the 
essential features of the short-period or Dutch-roll mode of the lateral 
motion were measured as a function of time and Mach number. 


The following quantities were then determined from the measurenents; 
period, damping factor, undamped natuiral circular frequency, lateral -force 
coefficient, time lag between sideslip angle and rolling angular accel- 
eration, and the an^jlitude ratios of rolling angular acceleration tc side- 
slip angle and lateral force to sideslip angle. An atten^)! was made to 
determine the amplitude ratio of yawing angular acceleration to sideslip 
angle, but this was not possible because of the erratic readings of the 
transverse accelerometer in the nose of the model. These erratic 
readings were due to the model experiencing both random and reguilar high- 
frequency oscillations superimposed upon the short -period osclllaticns in 
such a manner that the readings of this instrument were unreliable. 


The solution of the lateral equations of motion (eqs. 9 ho ll) by 
use of the time -vector method to determine the lateral stability deriv- 
atives is illustrated in figure 7 Doi’ hhe lateral oscillations following 
a rudder inpulse at M « 1.10. The lateral-force equation was solved 
first (see fig. 7(a)) to obtain the yawing angular velocity. The dis- 

ipI 

placement vector amplitude -|^ was plotted to the right with magnitude 

equal to unity in order to determine the scale of the vector polygon. 

The phase angle was considered to be 0°. The first vector plotted was 

mV ill vector where the amplitude ratio iiL 
‘IS 131 13| 

M 

quency (u)^ + aS)V^ , and the phase angle was determined from the 


the 


was determined by- 


multiplying the unit vector 


damping angle ^ 90 ° + tan”^ 


by the undanped natural circular fre- 
PP 

See eq-uation 4. This vector was plotted 
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in the positive direction of p. To this vector was added the 


|cyI IpI 
|Pl ' I PI 


vector and, since - j^ - is negative, it is plotted in the positive p 
direction. The vector - — a M 

was added next, where the amplitude 

IS IPI 


M 


ratio dJtll vas determined from the measured amplitude ratio 

ipi 




PI 


by 


j£ 


dividing py the undamped natiiral circular frequency, and the phase 

angle ^ was determined from the dan 5 >ing angle ^90° + tan“^ Note 


for this configuration since a is negative the vector a-L^ is 

<lS IPI 

plotted in the positive direction of The vector length now required 

I * I 

to close the vector polygon is the anplitude ratio — -pLL and the phase 

as IPI 

angle 0 ^^ is the angle measured between the \|r vector and the p vec- 
tor. Inasmuch as the yawing angular velocity was known, the yawing 
angular acceleration was determined by use of the basic properties of 
time vectors. 


After these quantities have been determined, either the rolling- 
moment or yawing-moment equation may be solved in a manner similar to 
the lateral -force equation. In both of these equations, the inertia 
vectors are known completely, whereas the aerodynamic vectors are known 
only in direction, so that two of these aerodynamic vectors must be 
solved for as a function of the third. Generally, it is best to esti- 
mate the arrplitude of the smallest vector so that greater acciiracy is 
allowed and then determine the other two by closing the vector polygon. 
However, transonic wind-tunnel values of C 7 are available for this 

'■P 

configuration from reference 9 j thus, values of and Cj were 

determined from the vector solution of the rolling -moment equation. 

(See fig. 7(ti)*) Values of ^*^n^ " ^ plotted in ijr direction, 

and were determined frcm the vector solution of the yawingHnoment 

equation as fxanctions of C 


“P* 

ass-umed to be -0.10, 0, and 0.10, and values of C- 


were 


(see fig. T(c).) Values of 

bp (s - ‘=np) 
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were deijermined for each of "these "values of order "to show tha"t 

C„ and (c^ - are linearly dependent upon C„ in the vector 

plots. It is seen from figure 7 (c) that Cn^ has little effect on Cn^ 

whereas it greatly affects The variation of ^n^ with Mach 

number shown later in the paper is presented for = 0 , whereas the 

variation of - C„. with Mach nmher is presented for all three 

n^ np 

values of C_ . 


RESULTS AND DISCUSSION 
Time History 


A portion of the time history of the lateral motions experienced hy 
the model due to a yaw distixrhance caused by an impulsive deflectior. of 
the rudder is shown in figure 8 where the variation with time of the 
rudder deflection, rolling angular acceleration, lateral-force coefil- 
cient, angle of sideslip, and angle of attack are presented. The appear- 
ance of the oscillations is that of a lightly damped sinusoid. 

Althoijgh the model was distvirbed essentially in yaw, an induced, 
pitching motion was recorded. It is believed that this angle -of -at"tack 
change is small enough to have a negligible effect on the lateral mctions 
The possible effects of inertia coupling experienced by some rocket- 
propelled models (see, for example, ref. lo) have been checked for "this 
conf ig"uration and the effects were found to be small. 

From the time history of the lateral motions of the model, the 
general characteristics of the Dutch-roll oscillation were determined. 
These characteristics are sho"wn in figure 9 where the period, total 
damping factor, and undamped nat"ural circular frecpuency are presented. 
Note that the period and the undamped natural circiilar frequency show 
the us"ual variation with Mach n"umber, that is, decreasing and incres.sing, 
respectively, with increasing Mach n"umber whereas the total damping fac- 
tor is low throughout the Mach n\xmber range of the tests and near 
M = 0.95 it is unstable. The period and total damping factor were deter 
mined from the oscillations in sideslip angle, lateral— force coefficient, 
and rolling angular acceleration, and the undamped natural circvilar fre- 
quency was determined from the faired variation of the period and total 
damping factor with Mach niimber. 
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Trim Characteristics 

The variations of the model trim characteristics with Mach nmber 
shown in figure 10 as trim sideslip angle, trim lateral-force coef- 
ficient, and trim angle of attack. There were no abrupt trim changes; 
however, the model changed trim somewhat throughout the test Mach n\mi- 
ber range. The trim change experienced in angle of attack was similar 
to the trim angle -of -attack variation reported in reference 1. The 
trim normal-force coefficient varied from 0.04 at M = 0.88 to O.O 5 at 
M = 1.1. 


Ang^litude Ratio and Phase of Rolling Angular Acceleration 

to Angle of Sideslip 

The variation with Mach number of the amplitude ratio and phase of 
rolling angular acceleration to angle of sideslip due to a yaw disturb- 
ance is shown in figure 11. It is important to note that these dynamic 
characteristics are for the model and the f-ull-scale airplane may or 
may not have the same phase or amplitude ratio since the mass and iner- 
tia characteristics of an airplane and a model would be considerably 
different . 


Sideslip Derivatives 

Cross plots .- A typical cross plot showing the variation of lateral- 
force coefficient with angle of sideslip is shown in figure 12. Note 
that, although the model was symmetrical and the rudder did not float 
when undeflected, zero lateral force did not occ\ir at zero angle of 
sideslip. This could be possibly due to an instrument shift during the 
flight. Since there was very little scatter in the data points from the 
faired variation of Cy against 3 and very little hysteresis thro\jgh- 

out the Mach mmiber range of the test, it is felt that the sum of the 
derivatives Cy and Cy is zero. 


Lateral -force derivative .- From the cross plots of Cy against p 

the lateral-force derivative Cy was determined and the variation of 

P 

Cy with Mach number is shown in figure 15(a) . This variation is similar 

p 

to the variation of the slope of the lift curve for wings having moderate 
aspect ratio, low sweep, and thick airfoil sections; thus, the main con- 
tribution to Cyp is the vertical tail. Furthermore, it is believed 

that the high value of Cy^ near M = 1.0 is greatly due to the induced 
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sidewash across the vertical tail. Tail-on and tail-off data at low 
and high speeds indicate a favorable induced sidewash effect and it is 
believed that the high value of Cy near M = 1.0 measured in the 

present test could be due to this favorable sidewash effect but of 
larger magnitude. 

Directional stability derivative .- The variation of the directional 
stability derivative C 


■^nr 


of the vector analysis is shown in figure 15(b) 
for 


with Mach number as determined from the results 

The data are presented 


an increase in 


with increasing Mach number. Calculation of the 


C„ =0 as explained in the analysis section and indicate generally 
T 

s 

isolated tail contribution to C„ based on the Cy data of fig- 

"P ^P 

ure 13(a) indicates that the isolated tail contribution is considerably 
greater than the data of figure 15(b) throughout the Mach number range 
of this test. It is believed then that the wing -body contribution tc 
C is destabilizing and of rather large magnitude. Also shown in fig- 

ure 13(b) is a conparison of C^j as determined from the s ingle -degree - 

of -freedom method. (See ref. ll) . The single -degree -of -freedom method 


neglects the product -of -inertia term and assumes 




in the yawing- 


moment eqixation; thus, Cj^^ is determined directly from the period of 
the yawing oscillation. Note that the agreement between Cj^ as deter- 
mined from the vector analysis and the single-degree-of -freedom method 
is good. This result is due to assimiing C^^ =0 in the vector analysis 

and also due to the fact that for this configuration the product of iner- 
tia is very small. It appears then that for configurations for which the 
product of inertia is very small and for which C„ is nearly equal to 

T 

zero, the single-degree-of -freedom method allows the determination of 
Cnp values that are good approximations to the actual values of Cj^ . 

These data then f\irther substantiate the results obtained in reference 12 . 
The C„ data of the present test determined from the s ingle -degree-of- 
P 

freedom method and the data of reference 1 do not show good agreement 
when these data are converted to the same center -of -gravity location. 

This resiilt is primarily due to two reasons. The data of reference 1 
were determined from rather random oscillations in yaw that were induced 
from pitching oscillations and also the moment of inertia in yaw of ttie 
model of reference 1 was estimated, not measured. 
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Effective dihedral derivative.- As discussed in the analysis section. 


values of the effective dihedral derivative C^ 


were determined from 


'P 


vector analysis by knowing values of the damping -in -roll derivative . 

P 

The variation of C7 with Mach number is shown in figure 13 (c) for the 

low lift coefficients of this test. For con^leteness, the data of 

reference 9 are included as figure l 4 . It is felt that the values of 

shown in figure 13(c) are primarily due to the vertical tail and that 
wing-body interference effects almost cancel the contribution of the wing 


dihedral to Ci 


'P 


Calculations of the isolated vertical-tail contribution 


to Cv tend to confirm this effect. Also high-speed wind-tvinnel data 

‘p 

indicate this effect; however, low-speed wind-tunnel data do not. 


Moment Derivatives Due to Yawing 


Rolling -moment-due-to-yawing derivative . - The variation of the rate 

of change of rolling -moment coefficient with yawing -angular -velocity 

factor C^ with Mach mmiber is shown in figure 15(a) . Values of C7 
tr I'r 

were determined from the vector analysis by knowing values of C7 , and 

data indicate that for the low lift coefficients of this test C7 is 

^r 

negative at M = 0.9 and increases positively with increasing Mach num- 
ber to a large positive value at M = 1 . 1 . 


Damping -in-yaw derivative .- The variation of the damping -in-yaw 
derivative C„ - C„, with Mach nimiber is shown in figure 15(b) where 

117. Xlp 

values of C„ - C„* are plotted for three values of C„ . These C„ 

n^ Up ^ np Up 

values were chosen in such a way that they wotild cover a range of possible 
values for the configuration. As may be seen, Cj^^ - C^j* is linearly 


dependent upon C„ ; thus, from the vector analysis, 
T 


obtained for any nimiber of assumed values for C 




Cn^ - Cn* may be 
The data plotted in 


figure 15(b) indicate that 0 ^^ - Cjj* is positive (vinstable) over a 


small region of Mach number near M = 0 . 95 j thus, possible dynamic insta- 
bility of the Dutch-roll oscillation was indicated. Actually as shown 
in figiire 9(b) , dynamic instability was recorded in the time histories 
of the lateral motions over the same Mach number range. This instability 
may be due to separated flow on the thick vertical tail at these Mach 
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numbers. A comparison of the data in figin’e 15('b) with the total- 
dan5)ing-factor data in figiire 9 indicates that for this configuration 
is positive. 


Another interesting point can be seen from figxrre L5(‘b) . Since 
the total damg)ing of the Dutch-roll oscillation is constant at any 
particular Mach nimber, increasing positively resiilts in values: 

of - C^. obtained from the vector diagrams that are less negatJ.ve 

Bj. np 

whereas decreasing C„ resiilts in values of - C^. that are more 


"P 


P 


negative so that the beneficial contribution of positive values of C 


“P 


to the dan^jing of the Dutch-roll oscillation becomes rather apparent, 
This effect is also shown in the vector plot of figure T(c). 


Con^jarlson of Present Data With Other Data and With Estimates 

Tests.- A comparison of the sideslip derivatives as obtained from 
the present test and the wind-tunnel tests of references I3 to 16 are 
shown in figure 16. Although there are rather large regions of Mach 
number where the sideslip-derivative data are not available, it is felt 
that the data from the present test together with the wind-tunnel test 
data may be \ised to determine faired variations of the sideslip deriA'a- 
tives with Mach number over the range of M = 0,l6 to M = 2.52. U. 
is also believed that the vector -analysis technique allows the deter- 
mination of the sideslip derivatives to the same order of accuracy as. 
does the wlnd-ttinnel technique since data from two wind tunnels at 
almost the same Mach number show at least a 10 -percent disagreement. 


Estimates.- A comparison of the lateral-stability -derivative dai.a 
of the present test and the estimated lateral stability derivatives i’or 
this configuration (ref. I7) is shown in figures I7 and I8. 


The coD5)arison of the sideslip derivatives is shown in figure I7. 
The differences between the values of Cy may be attributed to much 

larger values of sidewash than were estimated whereas the differences 

between the values of C, may be attributed to a stronger wing- 

ip 


fiiselage interference effect than estimated. Estimated values of 

taken from reference I7 were determined primarily from the test of ref- 
erence 1 where C„ was determined from the single-degree -of -freedom 
np 

method. The agreement between these data is good below M = 0.95 wHle 
the agreement at the higher Mach numbers is poor. As mentioned in a 
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previous section of this paper ^ the data of reference 1 were determined 
from random oscillations and an estimated moment of inertia in yaw so 
that the poor agreement in the C^^ data may he due to this fact. 


A comparison of the moment derivatives due to yawing is shown in 
figinre I8 where the agreement is poor. The C^^, data that were esti- 
mated do not show the abrupt variation with Mach number that the data 
of the present test indicate. The estimated values of Cj^ are for 

a = 0 whereas the data for the present test are for angles of attack 
that vary from a = - 1 . 2 ° at M = 0.9 to a = - 1 . 6 ° at M = l.lj how- 
ever, it is felt that the poor agreement is not due to this difference 
in angle of attack, although it is possibly due to the method used to 
determine the estimated values. The differences in the danping-in-yaw 
data may be attributed to the fact that the Cn^ contribution was not 

estimated in reference I7. 


Model High-Frequency and Random Oscillations 

Thro\;ighout the flight the instrimient telemeter system recorded 
high-frequency oscillations that were superln^josed on the low-frequency 
Dutch-roll oscillations. A portion of the telemeter record showing 
these oscillations is presented in figure I9. The high-frequency oscil- 
lations were rather regular in natvire from the peak Mach number of 
M = 1.18 to a Mach number of about M = O.96 and had a frequency of 
about 90 cycles per second throvighout this Mach number range. As stated 
in the preflight test section of this paper, shake tests were performed 
on the model. From these tests the first bending mode of the horizontal 
tall was determined to have a frequency of 86 cycles per second j however, 
the frequency of the vertical tail was not measured. The model of ref- 
erence 1 did not experience these high-frequency oscillations and, since 
the only difference between the model of reference 1 and the model of 
the present test was in the construction of the tail sections (weaker 
for present test) , it is believed that the model of the present test 
experienced flutter on some part of the empennage in the Mach number 
range from M = I.I8 to M = O.96. 

Below M = 0.96 the oscillations were random in nature j this condi- 
tion is believed to indicate that the model experienced a form of buffet 
which may be explained as follows: The horizontal tail was mounted low 

on the vertical tail near the section of the fuselage incorporating a 
large boattail so that the region of the model bounded by the vertical 
tail, lower surface of the horizontal tail, and the f-uselage represents 
a region where the static pressure may be considerably less than the 
free-stream static pressure and the flow may become unsteady. The 
results of reference 18 indicate that the static pressure in the region 
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discxissed above should be about O .5 of the free-stream static pressure 
from M =1.0 to M = 0.9 and that a low-lift buffet is possible in 
the region of M = 0.99 to O .87 because of the interference of the hor- 
izontal tail and body. It is believed then that the model experienced 
a low -lift buffet in the Mach number region below M = 0.95* 


CONCLUSIONS 


From the results of the flight test to determine the lateral sta- 
bility characteristics of a rocket-propelled model of a supersonic elr- 
plane configuration employing a sweptback wing having circular -grc 
sections, the following conclusions are indicated: 

1. The time-vector method applied to the recorded Dutch-roll tian- 
sient oscillations provided a useful method for the determination of 
the lateral stability derivatives. The results as obtained from the 
vector analysis indicate that, as the Mach number increases from O .89 
to 1 . 10 , the directional stability and the effective dihedral incre£.sed, 
the Hamping in yaw was low and over a small region of Mach number we.s 
unstable, and the rate of change of rolling -moment coefficient with 
yawing-angular -velocity factor was negative at a Mach number of O .89 
^nii increased positively to a large positive value at a Mach number 
of 1 . 10 . 


2. The lateral -force derivative increased with increasing Mach num- 
ber up to a Mach number of O .98 and then decreased to the limit Mac]:, 
number of the test. 

5 . For this configuration the lateral -force derivative was rather 
large, whereas the directional stability derivative was lowj thus, £. 
powerful sidewash effect on the vertical tail and a large wing-body 
destabilizing effect were indicated. 


k. The model experienced an empennage flutter in the Mach numter 
region above a Mach number of O .96 and a low-lift buffet below M = O. 96 . 


Langley Aeronautical Laboratory , 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January I 7 , 195^. 
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TABLE I.- GEOMETRIC Al® MASS CHARACTERISTICS OF THE MODEL 


Wing; 

Total included area, sq. ft 5 •5° 

Aspect ratio ^ 

Sweep of quarter-chord line, deg ^ 

Taper ratio 

Mean aerodynamic chord, ft 1.22 

Airfoil sections normal to quarter-chord 

line (max. thickness at 0.50c) 10 percent circular arc 

Horizontal tail: 

Total included area, sqft 0.958 

Aspect ratio 5*72 

Sweep of quarter -chord line, deg hO 

Taper ratio 0.5 

Airfoil section parallel to fiiselage reference line . . MCA 65-OO8 

Vertical tail: 

Area (exposed), sq ft O.825 

Aspect ratio (based on exposed area and span) I.16 

Sweepback of quarter chord, deg 55 

Taper ratio 0.557 

Mean aerodynamic chord, ft 1.02 

Airfoil section, root MCA 27-010 

Airfoil section, tip MCA 27-008 

Fuselage; 

Fineness ratio (neglecting canopy and fairings) 9*58 

Mass: 

Model weight, lb 190.5 

Moment of inertia in yaw, I2, slug-ft^ l4.5 

Moment of inertia in roll, Ij^, slug-ft^ 4.0 

Inclination of principal axes, e, deg (down at nose) .... 4 

Product of inertia, Ixz,^ slug-ft^ 0.752 

Center -of -gravity position, percent c' 20 

Wing loading, Ib/sq ft 54.2 
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TABLE II.- ESTIMATED ACCURACY OF BASIC QUANTITIES 

I 

' [Values shown are positive or negative quantities] 







Acciiracy of 

- 






Mach 

w, 

percent 

I7 

percent 

percent 

(a) 

M, 

percent 

percent 

percent 

M 

'^T/g 

deg 

a, 

deg 

Q 

number 

percent 

ipr 

percent 

IPI ' 

percent 

sec“^ 

1.1 

0.5 

2.0 

4.0 

8.0 

1.0 

2.0 

2.6 

5.0 

2.0 

4 

0.5 

0.1 

0.9 

0.5 

2.0 

4.0 

8.0 

1.7 

5.5 

2.6 

3.0 

2.0 

4 

0.5 

0.1 


®'This error is due primarily to estimated accuracy of inclination of principal axis, 


J 


I 
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Figure 1,- Sketch showing stability- and body-axes system. Each view 
presents a plane of the axes system as viewed along the third axes. 
Arrows indicate positive direction of forces^ moments, and angles. 
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Figiire 2.- Ang\ilar relationships in flight. Stability- and body-axes 
system. Arrows indicate positive direction of angles. t] = a - e. 
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(a) General arrangement of supersonic airplane configuration. 

Figure 3*- Drawings showing dimensions and characteristics of supersonic 
airplane configuration. All dimensions are in inches unless noted. 
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(b) Details of wing and horizontal stabilizer. 


Figure Continued 
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(a) Three-quarter rear view. 


L-81369.1 



(b) Three-quarter front view. 


L-81568.1 


Figure !<■.- Photographs of model of supersonic airplane configiaration. 


CCNFIDENTIAL 



NACA RM L55A31 


CONFIDENTIAL 


33 



(c) Model-booster combination on launcher. L-82i|88«l 
Figure 4.- Concluded. 
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Figure 6.- Variation of Reynolds nuniber with Mach number. 
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(b ) Rolling -moment equation . 


Figure 7 .- Continued. 
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Figure J.- Concluded. 
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- Portion of time history of flight. M = 1.10. 
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(b) Total damping factor. 



(c) Undamped natural clrciilar frequency. 

Figure 9.- Variation of period of lateral oscillations, total damping 
factor, and undamped natural circular frequency with Mach nimiber. 
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Figure 10.- Variation of trim characteristics with Mach number. 
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Figure 11.- Variation of amplitude ratio -j-^ and phase angle with 

|P I 

Mach number. 


m 

TFT 


CONFIDENTIAL 



NACA RM L55A31 


CONFIDENTIAL 


43 



P,deg 


Figure 12.- Typical variation of lateral-force coefficient with sideslip 

angle at M = 1.10. 
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(a) Lateral-force derivative. 



(b) Directional stability derivative. 



(c) Effective dihedral derivative. 

Figure 15*- Variation of sideslip derivatives with Mach nionber. 
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Figure l4.- Variation of damping-in-roll derivative with Mach number as 

obtained from reference 9 . 
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(a) Rolling -moment-due -to-yawing derivative. 



(b) Damping -in-yaw derivative. 

Figure 15*- Variation of moment derivatives due to yawing with Mach number. 
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Figure 16 .- Comparison of sideslip stability derivatives from various tests. 
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Figure 17.- Comparison between experimental and predicted variation cf 

sideslip stability derivatives. 
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Figure l8. - Comparison between experimental and predicted variation of* 

moment derivatives due to yawing. 
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Figure I 9 .- Portion of telemeter record showing high-frequency and randDm 

oscillations. 
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